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I In this paper, the degenerate ground states of Z2 topological order on a plane with 

O 

O ' holes (the so-called surface codes) are used as the protected code subspace to build 
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a topological quantum computer by tuning their quantum tunneling effect. Using a 
designer Hamiltonian - the Kitaev toric-code model as an example, we study quantum 
QQ I tunneling effects of the surface codes and obtain its effective theory. Finally, we show 

how to do topological quantum computation including the initialization, the unitary 
transformation and the measurement. 
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(N ■ I. INTRODUCTION 
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^ ' Quantum computers are predicted to utilize quantum states to process tasks far faster 

than those of conventional classical computers. Various designs have been proposed to build 
^ , a quantum computer, such as manipulating electrons in a quantum dot or phonon in ion 



traps, cavity QED, nuclear spin by NMR techniques. DiVincenzo pointed out that in order 



to build a real quantum computer[l], the following five criterions should be satisfied : 1) 
Scalability of extendible qubits; 2) Initialization (Creation of highly entangled states); 3) 
Local operations on multi-qubit; 4) Measurement of entangled states; 5) Long decoherence 
time. In particular, the fifth criteria (low decoherence condition) becomes a trouble to reach 
the goal. 

Recently, people find that it may be possible to incorporate intrinsic fault tolerance into 
a quantum computer - topological quantum computation (TQC) which has the debilitating 
effects of decoherence and free from errors. The key point is to store and manipulate 
quantum information in a "non-local" way, namely, the "non-local" properties of a quantum 
system remain unchanged under local operations. An interesting idea to realize fault-tolerant 
quantum computation is anyon-braiding, proposed by Kitaev|2|, |3||. He pointed out that the 
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degenerate ground states of a topological order make up a protected code subspace (the 
topological qubit) free from error 0, ^. 



Topologica 
paradigm [£ 



order is a new type of quantum orders beyond Landau's symmetry breaking 
Q]. People know that there are two types of topological orders in two di- 
mensional S = 1/2 spin models - non-Abelian topological ordered state and Z2 topological 
ordered state. Those topological ordered states may appear in frustrated spin systems or 



dimer models 



0, Q, 



Because of full gapped excitations, these topological orders are 



robust against any perturbations, even those perturbations that break all symmetries. 

In non-Abelian topological orders, the elementary excitation are non-Abelian anyon 
with nontrivial statistics. Now people focus on realizing TQC by braiding non-Abelian 
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anyons|3l. 1141]. The degenerate states undergoes a nontrivial unitary transformation when a 
non-Abelian anyon moves around the other. One can initial, manipulate and measure the 
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degenerate ground states wit 



issue recently [15 
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ian anyonsjsl, [ij], which has become a hot 
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On the contrary, Z2 topological order is the simplest topological ordered state with three 
types of quasi-particles: Z2 charge, Z2 vortex, and fermions IJ] . Z2 charge and Z2 vortex 
are all bosons with mutual vr statistics between them. The fermions can be regarded as 
bound states of a Z2 charge and a Z2 vortex. In the last decade, several exactly solvable 
spin models with Z2 topological orders were found, such as the Kitaev toric-code model 2], 

on a honeycomb lattice [3|]. In 
Ref. 261], an alternative way to design TQC is proposed by manipulating the toric codes - 
the protected code subspace of Z2 topological orders. To manipulate the degenerate ground 
states we can tune tunneling by controlling external field on spin models. However, to build 
a real quantum computer, multi-qubit is necessary. It is indeed a challenge to realize a spin 
model on a manifold with higher genus in experiments. To solve this problem, in this paper 
we design TQC by manipulating the surface code (See detail in the main content) rather 
than toric code, since it is more easy to make a hole in a surface. 

The paper is organized as follows. In Sec. II, we study the properties of the surface code. 
In Sec. Ill, an effective theory of the surface code in Z2 topological orders is formalized. 
By using the Kitaev toric-code model as an example, we demonstrate how to control the 



surface code by tuning the tunneling of the degenerate ground states 



In Sec. 



IV, based on the effective theory, the TQC is shown including the initialization, the unitary 



FIG. 1: The scheme of the Kitaev toric-code model 



transformation and the measurement. Finally, we draw the conclusions in Sec. V. 



II. DEGENERATE GROUND STATES AS SURFACE CODES 



In the paper, we demonstrate the TQC in topological order by using a designed model 
the Kitaev toric-code model (an effective model of the Kitaev model on two dimensional 
hexagonal lattice) as aji example [2], s], J]. Here the Hamiltonian of the Kitaev toric-code 



model is described by 



where 



(1) 



I Seven 



is odd 



with > 0. cTj^'^'^ are Pauli matrices on sites i. See the scheme in Fig.l. 

Firstly we study the ground state degeneracy Q of the Kitaev toric-code model. The 
ground state is a Z2 topological state that is denoted by = = +1 at each site with 
energy, 

Eo = -gM (2) 
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271 ]. Under the periodic boundary condition 



(on a torus), the degeneracy Q is dependent on A/" : Q = 4 on even- by-even (e * e) lattice, 
Q = 2 on other cases (even-by-odd (e * o), odd- by-even (o * e) and odd- by-odd (o * o) 
lattices) 0, 12, 13, 27]. In addition, on a manifold with high genus (x > 1), Q becomes 4?^ 



on e * e lattice and 4-*^ — 2 on other cases. Kitaev have noted the degenerate ground states 
of the Kitaev toric-code model on a torus as the toric code. On the other hand, for the 
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FIG. 2: Mapping the qubit of a hole onto a pseudo-spin 

ground states of the Kitaev toric-code model on a surface with open boundary condition, 
the ground state degeneracy now becomes 

Q = 2" (3) 



(without considering the edge states) |28l.l29l|. Here n is the number of holes. In the following 
part we call the degenerate ground states of the Kitaev toric-code model on a surface with 
open boundary condition the surface codes. 

To classify the degeneracy of the ground states (the surface codes) , we define three types 
of closed string operators Wc{C), W,,{C) and Wf{C) = Wc{C)W^{C). Here Wc{C) = He < 
(or Wy{C) = Ylc^i) products of spin operators along a loop C connecting even- 

plaquettes (or odd-plaquettes) of neighboring links, with C denoting closed loops. It is 
obvious that the three types of closed string operators {Wc{C), Wv{C), Wf{C)) correspond 
to three quasi-particles {Z2 charge, Z2 vortex, and fermions) respectively. One can easily 
check the commutation relations between the closed string operators and the Hamiltonian 

[H, W,{C)] = [H, W,{C)] = [H, Wf{C)] = 0. (4) 

In particular, for the surface codes, we can define two types of special closed string 
operators, Wy{CA) and W/(Cb). Here Ca denotes a closed loop around a hole (labeled by 
an index I) and Cb denotes a loop from the hole to the boundary of the system. Although 
Wf{CB) is not the original closed string operator, its topological properties are the same. 
Due to the anti-commutation relation between Wy{CA) and W/(Cb), 

{Wy{CA),Wf{CB)} = 0, (5) 
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FIG. 3: Mapping the qubits of a line of holes onto a pseudo-spin chain 

we may identify Wu(Ca) and Wf{CB) as pseudo-spin {S = |) operators rf and rf , respec- 
tively. The ground states become the eigenstates of r^^. Then one has two degenerate ground 
states (denoted by | mi)) for the case with a single hole. For mi = 0, we have 

Ti I m;) =1 mi), (6) 

and for mi = 1 we have 

\mi) = - \ mi). (7) 

Physically, the topological degeneracy arises from presence or the absence of tt flux of Z2 
vortex through the hole (See Fig. 2). The values of m; reflect the presence (m; = 1) or the 
absence {mi = 0) of the tt flux in the hole. 

Furthermore for the Z2 topological order with n holes, the degenerate ground states have 
2" fold degeneracy as Eq.(3). The 2" degenerate ground states can be mapped onto a 2"-level 
quantum system of pseudo-spins fi by the following correspondence, 

|m; = 0)^|T),, \mi = l)^\i)i. 

We now get 2" degenerate ground states of a chain of n holes. The 2" degenerate ground 
states are denoted by 

I mi, ...,mi...,mn) =\ mi) (g) | mi)...® \ mn), 

where the two degenerate ground states of each hole are denoted by | m^), m; = 0, 1. In the 
following parts we use the degenerate ground states on a surface with a chain of holes to do 
TQC (See Fig.3). 
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III. PERTURBATIVE APPROACH OF QUASI-PARTICLES 

In this section we study the properties of the quasi-particles. In this solvable model, Z2 
vortex is defined as = — 1 at even sub-plaquette and Z2 charge is Xj = — 1 at odd sub- 
plaquette. The mass gap of Z2 charge and Z2 vortex is 1g. In particular, the bound states 
of a Z2 charge and a Z^ vortex on two neighbor plaquettes obey fermionic statistic. All 
quasi-particles in this model have fiat bands. The energy spectra are = = 2g for Z2 
vortex and Z2 charge, Ef = Ag for fermions, respectively. In other words, the quasi-particles 
cannot move at all. 

Under the perturbation 

i ill 

the quasi-particles begin to hop. The term ^ drives the Z2 vortex without affecting 

i 

fermion and Z2 charge. For a Z2 vortex at i plaquette Xi = —1, when af acts on z + site, 
it hops to i + ix — iy plaquette denoted by Xj+e^_ey = — 1, 

X= +1 ^ Xj = —1, Xj_|_e^_ej^= +1 Xj_|_e^_ey = —1. 

Moreover, a pair of Z2 vortices at i and i + ix — Cy plaquettes can be created by the operation 
of af , 

X= +1 — ^ Xj = —1, Xj_i_g^_gj^= +1 — ^ Xj_i_g^_gj^ = —1. 

Similarly, the term ^ drives the Z2 charge without affecting fermion and Z2 vortex. 

i 

In particular, there exist two types of fermions : the fermions on the vertical links and the 
fermions on the parallel links. The term ^ cxf drives fermions hopping without affecting 

i 

Z2 vortex and Z2 charge : the fermions on the vertical links move vertically and the fermions 
on the parallel links move parallelly. That means both types of fermions cannot turn round 
any more. 

To describe the dynamics of the quasi-particles we use the perturbative approach in 



Ref. 
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20 



21 
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26| . In the perturbative approach, the spin operators are represented 
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by hopping terms of quasi-particles, 



Here 0)^^^ ( a = 1,2 ) are the generation operator of vortex, Z2 charge and ^^ia ( 
a = 3, 4 ) are the generation operator of fermions, respectively. i\ denotes the position on 
even sub-plaquette and ii denotes the position on odd sub-plaquette. ^3 denotes the position 
on the vertical links and 24 denotes the position on the parallel links. In addition, one should 
add a single occupation constraint (hard-core constraint) as 



^) = 1^) or 0^1^) I ^) (10) 



where denotes quantum state of the Kitaev toric-code model. Therefore, by the per- 
turbation method, the Hamiltonian H = Hq + H' can be represented by generation (or 
annihilation) operators of quasi-particles. 



+4^? ^Us^^^^^ + 5^ <PhAM (11) 



14 



and 



H' - -h^ E i4,n<Pi,n + <PlA,n) (12) 
-h' Xl('^U'^2,i2+0U0L2) 

4,i3*^3,j3±e^) 



-^'5^('^3,i3'^3,to + 

*3 

Y{(j){i^(j)4M±ey + 0l,i4 0li4±e,) 
(*3J4> 

In the following parts, we consider only the perturbation as H' = h^^cr^ + h^Yl^i- 

i i 

The perturbative Hamiltonian of the quasi-particles becomes + Hc + Hj where 
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«i ihJi) 
He = 2^X]02,i2</'2,i2 

is 14 43 

5^(01,^4 </'4,i4±e, + 4,^AM±ey) + h.C. (13) 

u 

That means Z2 charge cannot move any more. 

For the perturbative Hamiltonian on Ly x square lattice {Ly^ are all even integers), 
we obtain the dispersion of quasi-particles. The energy of Z2 vortex is given by 



e,^ = ^{^,^+2gy-e,, (14) 

where 

^fc, = 2h''[cos{k^ + ky) + cos{k^ - ky)]. (15) 

Here kx and ky are the wave vectors (The lattice constant has been set to be unit). The 

energy gaps of Z2 vortex is obtained as 2g^l — On the other hand, the energies of 
fermion is given by 

Sks = 4^7)2 + 4(/z2^ COS fc,.)2 (16) 



^fcs = \/{^k4+4^9y + 4^ihy cos kyY 



2hv 

a 

dispersion of Z2 vortex and fermion by tuning the external field. 



The energy gaps of fermion is obtained as Ag^jl — Thus one may manipulate the 



IV. EFFECTIVE PSEUDO-SPIN MODEL OF SURFACE CODES 

It is known that the degenerate ground states of Z2 topological orders have identically 
energy in thermodynamic limit. However, in a finite system, the degeneracy of the ground 
states is (partially) removed due to tunneling processes, of which a virtual quasi-particle 
moves around the holes before annihilated with the other . In general cases, one will 

get very large energy gaps (a = 1,2,3,4) for all quasi-particles and very tiny energy 
splitting of the degenerate ground states AE'. That is /S.E <ti nia- Based on this condition 
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FIG. 4: Tunneling process of Fermions from one hole to the boundary 
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FIG. 5: Tunneling process of Z2 vortex moving around one hole 



{AE <^ rria), we may ignore excited states with E > rria and consider only the topological 
degenerate ground states. Thus we get a 2" system as the effective pseudo-spin model of 
the surface code. In the foUowings, we will derive this model step by step. 

In face, the closed string operators Wy{CA) and IV/(Cb) can be considered as quantum 
tunneling processes of virtual quasi-particle moving along the loops. Let us take the quantum 
tunneling process of fermions as an example : at first a pair of the fermions is created. One 
fermion propagates around the hole driven by the operator af and then annihilates with the 
other. Then a closed string of af is left on the tunneling path behind the virtual fermion, 
that is just a closed string operator Wf{CB)- Such a process effectively adds the vr flux to 
one hole and changes by 1. 

Firstly we calculate the the effective pseudo-spin model of single qubit (labeled by l-th 
hole) from tunneling processes. When a virtual fermion propagates from the boundary of 



the l-th hole to the boundary of the system, the quantum state 




turns into 





(17) 
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Such process is shown in Fig. 4. By a degenerate perturbation approach j^, 0, one may 
obtain the energy sphtting Ai? of the two ground states as 

6E = J2 {niH'i-^H'Ylih (18) 

j ^0 — -no 

which becomes 

SE = 2-^^^-^^ (19) 

where is the length of the shortest path of a fermion from the boundary of the /-th hole 
to the boundary of the system. See detailed calculation in Ref . jsojl . 

On the other hand, considering a virtual Z2 vortex propagating around the hole (shown 

IT)^ 

m, \ I \\\\ 

(20) 



in Fig. 5), the quantum states I I turn into 



The corresponding energy difference e of the two ground states is 



(21) 



where is the length of the shortest path of a Z2 vortex around the hole. 

Thus the dynamics of such a two-level quantum system (a single qubit) can be described 
by a simple effective pseudo-spin Hamiltonian 

Heff = ^(iT)zai/ + u)Kri/) + |(ir)Krb-ii)/aio 

= }xlTt + h\Tl (22) 

with hf = ^ and hf = |. 

Secondly we calculate the effective exchange interaction between two qubits. For simplic- 
ity, we consider only the perturbation as H' = '^i + X] '^1- Then from Eq.(13), there 

i i 

exist two different tunneling processes : virtual fermion propagating from the boundary of 
/-th hole to the boundary of (/ + l)-th hole, virtual Z2 vortex propagating around the two 
holes, respectively. See Fig. 6 and Fig. 7. Let us calculate the ground state energy splitting 
from degenerate perturbation approach. 

When a virtual fermion propagating from the boundary of /-th hole to the boundary of 
(/ -|- l)-th hole, the quantum states turn into 
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FIG. 6: Tunneling process of fermions from one hole to another 
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FIG. 7: Tunneling process of Z2 vortex moving around two holes 






Therefore, by the mapping, the pseudo-spin operator of the tunneling process of fermion 
corresponds to rf ®r;^]^. See Fig. 6. The energy splitting 6E of the ground states is obtained 
as 

SE = Iff,. = -^-^^ r (23) 

where Lyy is the length of the shortest path from the boundary of /-th hole to the boundary 
of (/ + l)-th hole. 

Similarly, when a virtual Z2- vortex propagates around /-th and (/ + l)-th holes, the 
quantum state turns into 





The pseudo-spin operators of the tunneling process of Z2 vortex correspond to Tj^CSTj^i. The 
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tunneling amplitude is 

where L^x is the length of the shortest path round both /-th and (/ + l)-th holes. Thus under 
the perturbation H' = af + erf, the total effective Hamiltonian of the exchange 

i i 

interaction becomes 

Kff ^ E(^M+i^^r+i + JtUiT^rl^i)- (25) 
I 

Finally, for a chain of n-hole, the degenerate ground states can be mapped onto a model 
of a n-pseudo-spin chain. By ignoring the next nearest neighbor coupling terms, the effective 
model of the Kitaev Toric-code model under the perturbation H' = h^'Y^ erf + h^'Y^ erf is 

i i 

naturally an anisotropy Heisenberg model 

Kff ^ $^(J^r+irrrr+, + Ji%,t,\%,) + Y,Chfr^ + hfrn. (26) 

I I 

The effective Hamiltonian in Eq. fl26l) indicates that the Kitaev Toric-code model is an exam- 
ple of so-called topological order with controllable dispersion of quasi-particles. For example, 
if one adds the external field along y-direction only encircling two holes, H' = W Y '^f ? the 

i 

effective model is reduced into 

Kff ^ Jfjrfr^^ + /i^r + hlrl (27) 

So one may adjust each parameters in the effective Hamiltonian by controlling the local 
distribution of the external field along special direction. 



V. TOPOLOGICAL QUANTUM COMPUTATION WITH SURFACE CODES 

To design a topological quantum computer, one needs to do arbitrary unitary operations 
on the surface codes. Then by adding the specific perturbations to the Kitaev Toric-code 
model, H\ one can change different quasi-particles' hopping and then manipulate the surface 
codes by controlling tunneling splitting of degenerate ground states. In this part we show 
the initialization, the unitary transformation and the measurement. 
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A. Initialization 



Firstly we will show how to initialize the system into the quantum state |f , t)- The 
basic idea is to polarize all pseudo-spins by adding an effective field along x-direction and 
then removal it slowly. This process will occur according to the Hamiltonian 

H'=h^it)J2at (28) 

i 

where h^{t) = /iQe~*°/*. At the beginning, there is a finite external field h^{t) ^ h^^t ^ —oo. 
At the time t = 0, the external field disappears, h^{t) = 0. From Eqj26l we get the effective 
pseudo-spin Hamiltonian of the surface codes 

Kff ^ E -JiUi^trU, + hmri (29) 
I I 

where Jf^f+i = ^^l^g^lZ-^ /if(t) = ^^^*^ix_i ■ Here Hq is positive and is an even 
number. Then if the system evolves adiabatically and continuously from high temperature 
to the ground state, after a long time, the final state can be a pure state of the topological 
order ||, |) which becomes the initial state prepared for TQC. 



B. Unitary operations 



Secondly we discuss how to do an arbitrary unitary transformation on the surface code 17|, 



231 ] ■ The key point here is that the unitary operations can be achieved by controlling the 
external field along particular direction within fixed times. 

A general pseudo-spin rotation operator of l-th qubit is defined by 

Ui{e, if, 7) = e-^^^'e-^^^i'e-^'^^' (30) 

where 7 = hfAt^, 9 = hfAtg and ip = hfAt^. One can use external field along different 
directions encircling only l-th hole to do TQC : firstly applying the external field along 
y-direction at an interval Atg. The effective Hamiltonian becomes Tics = hlr^ . Then, we 
swerve the external field along x-direction at an interval At^p. The effective Hamiltonian 
becomes TicS = hfrj^. Finally, the external field along y-direction is added at an interval 
At^. The effective Hamiltonian becomes Tiefr = h^Ti^. 
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Using such method, one can reach certain quantum operations demanded by TQC and 
have the abihty to carry out gate operations onto the 'local' qubit (Z-th hole) at will, 

|^)/ = «|T)/ + /9e^*li)/ (31) 

with a, (3 > (a^ + = 1). An example is the | gate of a local qubit (/-th hole), of which 
we have the general pseudo-spin rotation operator as 

Ui{7= I, 9=0,^ = ^). 

One can firstly applying a global external field along y-direction at an interval At^ = 
Then, we swerve the global external field along x-direction at an interval At^ = Simi- 
larly, one may design the global Hadamard gate as a special pseudo-spin rotation operator 
on each qubit, 

TT( ^ a ^\ 

Thus, in principle, people are capable of to do arbitrary unitary transformation on the 
protected subspace by controlling the external field on given regions (for example, a close 
loop around one or more hole). 



C. Measurement 



Thirdly we discuss the measurement of an arbitrary quantum states of the surface codes. 
The central point is to measure the expected values of pseudo-spin operators by observing the 
quasi-particles' interferences from Aharonov-Bohm (AB) effect. 

To determine a, [3 and of quantum state of the /-hole, | = a |t)« + /^e*'^ we 
need to observe both fermion interference and Z2 vortex interference. Fig. 8 is a scheme to 
show the AB interference. 

Firstly we detect the value of (rf) = \i rf \ to determine a and l3 by AB effect 
from Z2 vortex-interference. To observe the AB interference, we add a small external field, 
/i^ — > and = 0. Now Z2 vortex begin to hop. There exist symmetrical paths from both 
sides of the hole. For example, 71 and 72 shown in Fig. 8 are two symmetrical paths. Then 
the symmetrical trajectories will contribute to the transition amplitude Tjj according to : 

^M = kZ}r + |V^Mf + 2e|V'Z5-V'5-| (32) 
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FIG. 8: The scheme of the interference around Z-th hole. 

where '07^ and ■j/'/i the wave functions of Z2 vortex of the two trajectories. For the ground 
state II)/, e is unit. However, for the ground state with vr-flux inside the hole ||)/, we have 
e = —1. Then we can distinguish these two cases. For two symmetrical paths ^/'^j — '^Tj~ 
we get a probability for |t)/ with Tjj = 4tj and a probability for with Tjj = 0. 

On the other hand, one can detect the value of (rf ) = |; t'i \ to determine the 
parameter by observing fermion interference. To observe the AB interference of fermion, 
we add a small external field, ^ 0. The wave function of fermion has a periodic boundary 
condition from the hole to the boundary of the system for the ground state \y)i = \] 
)i+ and an anti-periodic boundary condition for the ground state = ^ |T)z— 
Then an arbitrary state a |t)/ + /^e**^ \ is re- written into 

yi + a/3cos0e^'^' |T')^ + ^ ^ - a/3 cos 0/3e^<^" (33) 



where 

rh' = flTctan ( 

' P COS (j) + a 

and 



arctan(— ) (34) 



,,, , sm0 . . 

(p = arctan(— ). (35) 

p cos (j) — a 

For two symmetrical paths (one from the hole, the other not), we get a probability (| -|- 
aP COS0) for |t'); with Tij = itl and a probability | — aj3 coscp for with Tij = 0. As a 
result, we determine the parameters a, (3 and of an arbitrary state | vac) = a \'\)i+l3e^'^ ||);. 

The situation becomes more complex for 2-qubit. For 2-qubit with 2 holes (Z-th and 
/ -|- 1-th), a general entangled quantum state is given by 

«i I T)m® \])i + a2e"^' |T)m® IT)« 
+a^e'^^ I T);+i® |i); + a4e^'^Mi);+i®lT)z 
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Fig9 (a) 
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Fig9 (c) 



FIG. 9: The scheme of the mterference around l-th. hole in (a), the interference around / + 1-th 
hole in (b), and the interference around two holes in (c). Z2 vortex cannot move in shadow region. 

Here and 0j are all real number. Because of the constraint, Xlj'^? ~ there are totally 
2(2^ — 1) = 6 independent parameters. One can detect the value of (rf ) , {ti^i) , {t[ti\-^) 
to determine (z = 1,2,3,4) from the AB effect of Z2-vortex-interference by 2^ — 1 = 3 
times measurements. As shown in Fig. 9, there are 3 cases of the Z2 vortex-interference of 
different paths : 7" and 72 are two symmetrical paths around the hole l; 7J and 72 are 
two symmetrical paths around the hole / — 1; 7J and 72 are two symmetrical paths around 
both hole / and hole / — 1. To do the observations, we only add a small external field in the 
regions without shadow, and h'^ = (In the shadow regions, there is no external 

field, = = 0). So the Z2-vortex is guided moving around given holes. Similarly, one 
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can detect the value of (rf) , (t^+i) , (Tj^T^+i) determine </)j (z = 2,3,4) from the AB effect 
of fermion-interference by 3 times measurements. 

For n-qubit with n holes, a general entangled quantum state is given by 

5;]«ie^^MT)i® IT)^® lT)m (36) 

i 

Here and (pi are all real number. Because of the constraint, Xlj*^? — ^5 there are totally 
real 2(2" — 1) independent parameters by changing the paths of quasi-particles. So people 
needs to do 2(2" — 1) times measurement to determine the entangle state by changing 
different paths of AB interferences. One can determine ctj from "'^"■^^^ times measurements 
of Z2 vortex-interference (to detect (rf ) , (rfr;^^) , (Tir[j^^Ti_^^ (WiT^) ) and (pi from 
^^'^^^'^ times measurements of fermion-interference (to detect (rf ) , (Tfr^^^^)), respectively. 
So one cannot determine all the parameters by the quasi-particles interferences if n > 3. It 
is still an unsolved problem to measure a general entangled quantum state of n-qubit with 
n ( n > 3) holes. 



D. Errors 



Finally we discuss the errors and the constraint on our proposal. 

Errors mainly come from the thermal effect. At finite temperature, Z2 vortices are ex- 
cited, their moving around the holes leads to errors, as causes "thermal hopping" from one 
degenerate ground state to another. The decoherence time tdo has been roughly estimated 
by the time to stretch a pair of Z2 vortices over a distance equal to the average inter-particle 
separation (This is because the energy gap of Z2 vortex is smaller than that that of fermion). 
At low temperature, one can estimate that tdc is about t^e^sT where t^ = is the time 
scale for Z2 vortex moving the length of tunneling pathes, Lp 3l|]. v* is the average speed 
of Z2 vortex, which is estimated by v* ~ where Mea — {2h^)~^ is the effective mass 



of Z2 vortex. 

So there must exist a crossover temperature T* from thermal hopping to quantum tun- 
neling. Above T*, the decay rate of the quantum states is determined by process of thermal 
activation, which is governed by the Arrhenius law, 

to^e"^. (37) 
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FIG. 10: The crossover temperature for an example of tunneling process of Z2 vortex. The scale 
of temperature is — |^ 

Therefore at high temperature the errors proliferate and one cannot get reliable TQC. Below 
T*, quantum tunneling processes dominate, the rate of which goes as where B is about 



5 ~ max(Lp In — Lpln— ^). (38) 



Ignoring the prefactor and equating the exponents, one obtains 



(39) 



Fig. 10 shows the crossover temperature via for a quantum tunneling process of Z2 vortex. 
Thus if the temperature is kept far below T*,T <^T*, one may do unitary operations safely. 

However, because of the errors from the stochastic fields, we still get in trouble on storing 
quantum information by the surface codes. To store quantum information, all quantum 
tunneling processes need to be suppressed as low as possible, 0. That means the 

external fields should be removed, = = 0. Now the estimation of the crossover 
temperature in Eql39] is invalid. Without external fields, due to the diverge effective masses 
of Z2 vortex, M^s 00, one may store the quantum information for arbitrary long time, 
tdo ^ 00. Whereas, the stochastic noise fields leads to a finite decoherence time t^e- Thus to 
get a long-lived quantum information, both the temperature and the stochastic noise fields 



should be suppressed below a threshold. See detail in Ref. 32 1. 
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VI. SUMMARY 

In this paper, we find an alternative way towards designing a quantum computer that 
may be possible to incorporate intrinsic fault tolerance. Using the Kitaev toric-code model 
as an example, we obtain the effective pseudo-spin model that can be mapped onto the 
anisotropic Heisenberg model of a pseudo-spin chain in external field. Then one may tune 
the parameters of the effective pseudo-spin model by controlling tunneling processes of the 
surface codes by applying external field along special direction on lattices. In particular, 
five criterions to build a quantum computer are satisfied : 

1. Scalability of extendible qubits : The qubits is the so-called surface codes (two degen- 
erate ground states of topological orders on a plane with a hole). So the quantum 
computer is just a line of holes in the topological order. 

2. Initialization (Creation of highly entangled states) : We may polarize the pseudo-spins 
by adding an effective field along x-direction and then removing it slowly. 

3. Local operations on multi-qubit : The unitary operations can do by controlling the 
external field along particular direction within fixed times. 

4. Measurement of entangled states : We may measure the expected values of pseudo-spin 
operators by observing the quasi-particles' interferences from AB effect. 

5. Low decoherence : Below the crossover temperature T*, the decoherence processes will 
be controlled as low as possible. 

Finally we discuss the realization of the Kitaev toric-code model. Because it can be 
regarded as an effective model of the Kitaev model on a two dimensional hexagonal lattice, 
one may realize the Kitaev model firstly. The realization of the Kitaev model has been 



proposed in an optical lattice of cold atoms in Ref. 
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34| and in Josephson junction array 



of a superconductor in Ref. [35|. So it is possible to design a quantum computer in these 
system in the future. 

In general, to design a topological quantum computer via quantum tunneling effect, there 
are four necessary conditions : 1) Z2 topological order as the ground states; 2) controllable 
dispersion of quasi-particles; 3) space with nontrivial topological structure (manifold with 
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high genus or multi-hole) ; 4) low temperature and low noise. Thus considering these condi- 
tions, one may do TQC based on other models. 
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